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Abstract 



The quasistatic approximation and equation-of-motion decoupling for the electron Green's functions are applied to trace the effect 
of electronic dispersion and electron correlations on the ferromagnetism of two-dimensional itinerant-electron systems. It is found 
Q_ithat next-nearest-neighbor hopping t' is of crucial importance for ferromagnetism formation yielding the magnetic phase diagram 
& which is strongly asymmetric with respect to half-filling. At small f in the vicinity of half-filling the ferromagnetic phase region is 



in 



restricted by the spin-density wave instability, and far from half-filling by one-particle (spin-polaron) instability. At t' close to f/2 



ferromagnetism is stabilized at moderate Hubbard U due to substantial curvature of the Fermi surface which passes in the vicinity 
of the van Hove singularity points. The results obtained are of possible importance for high~r c compounds and layered ruthenates. 
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1. Introduction 

During last two decades, the explanation of magnetic proper- 
ties of itinerant-electron compounds attracts substantial interest 
in connection with physics of layered systems. The highly- 
correlated copper-oxide high-temperature superconductors, 
e. g., La2-. t Sr x Cu04, demonstrate under doping both commen- 
surate and incommensurate antiferromagnetism [1], and the 
phase separation (PS) related to magnetic phase transitions __|. 
Moderately correlated one- and bilayered strontium ruthenates, 
Sr2_ A La A Ru04, Sr2_ A Ca j: Ru04 and (Sri_ A Ca A )3Ru207, demon- 
strate the interplay of (commensurate) ferro- and (incommen- 
surate) long-wave magnetic fluctuations |__, and large Wil- 
son ratio |4] indicating proximity to ferromagnetic instability, 
which is achieved under doping. Explanation of drastically dif- 
ferent magnetic properties of these two compounds with iso- 
morphic crystal structure appears to be an intriguing challenge. 

From a theoretical point of view, the physics of the layered 
compounds can be qualitatively captured within the one-band 
Hubbard model on the square lattice, basing on Cu or Ru d n - 
derived band and accounting for nearest (f) and next-nearest 
neighbor hopping (f') integrals. ARPES and x-ray studies 
suggest the discrepancy between these compounds originating 
from substantial difference in their values of t' and Coulomb 
(Hubbard) interaction U. 

A large enough value of density of states (DOS) at the Fermi 
level (which is usually provided by van Hove singularity, VHS) 
leads to ferromagnetic ground state within the Stoner theory 
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(mean-field theory of ferromagnetism). However, this ap- 
proach neglects both Neel antiferromagnetism and incommen- 
surate magnetic ordering. Including these orderings into the 
mean-field theory __| allows to trace the stability of ground- 
state ferromagnetic order with respect to spin-density excita- 
tions depending on t' . It turns out that at small ratio t' jt the 
ferromagnetic ordering is stable only at non-realistically large 
values of U jt (even not too close to half-filling), and at smaller 
U/t or close to half-filling the Neel or spiral magnetic order is 
stable. 

At larger ?' the van Hove singularity (i) provides the stability 
of the ferromagnetic state with respect to spin-flip one-particle 
excitations provided that the Fermi level lies in the vicinity of 
VHS, (ii) is well separated from the chemical potential posi- 
tion at half-filling where the most preferable magnetic order- 
ings are characterized by large wave vector. It turns out that 
the competition of ferromagnetic and long-wave spin ordering 
in the vicinity of van Hove filling is a main limiting factor for 
ferromagnetic ordering. A priori it can be conjectured that the 
critical value of U which is needed to stabilize ferromagnetic 
order is much smaller than in the case of small t' . 

Another effect of finite ?' is strong dependence of PS region 
into the Neel antiferromagnetic and spiral (at moderate U) or 
ferromagnetic (at large U) phases on the charge sign of current 
carriers __. 

Above-discussed approximations completely miss correla- 
tion effects. The consideration of the case of large U enables 
one to see the problem of ferromagnetism stability in new light. 
Okabe [6] applied a variational principle which provides some 
interpolation scheme between the cases of moderate and strong 
electronic correlations. It was found that at small number of 
current carriers (holes) ferromagnetic ground state can be un- 
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stable with respect to spin-density wave formation, whereas 
at large hole number with respect to individual spin-flip. In 
Ref. 1 7] this problem was considered within a Gutzwiller-type 
method for t' + 0. It was found that in the case of large t' jt the 
ferromagnetic order acquires an additional stability for n > 1 (n 
is the electron concentration). 

During last decade some weak- and moderate coupling stud- 
ies taking into account the correlation effects provided some 
progress in understanding magnetic properties for non-zero 
t' . The spin-density wave instability of ferromagnetically or- 
dered ground state was considered for rather large t' (?' < t/2) 
within the quasistatic treatment of magnetic fluctuations (the 
quasistatic approximation) |8[] . It was found that critical values 
of stability of ferromagnetic ordering substantially exceed those 
determined in the Stoner theory. The interplay of magnetic and 
electronic properties and their effect on the possibility of fer- 
romagnetic instability was studied by the functional renormal- 
ization group (fRG) technique (see [9] and references therein). 
Main result of these studies is that small values (f jt < 0.2) fa- 
vor the antiferromagnetic instability at Van Hove filling, mod- 
erate ones (0.2 < f' jt < 0.35) the aPwave supeconducting in- 
stability, and large f jt > 0.35 correspond to ferromagnetic in- 
stability. 

In the present paper we investigate the ferromagnetism prob- 
lem for different f within the quasistatic approximation and 
many-electron approach. 

2. Instabilities of ferromagnetic state: t' jt far from 1/2 

We start from the Hubbard model 

H = ^ fkc^cw + U ^ c^c n c] iCi[ . ( 1 ) 

ktr i 

with the bare electronic spectrum e k = -2f(cos k x + cos k y ) + 
At' (cos k x cos k y + 1). 

To treat the spin- wave instability of the ferromagnetic ground 
state we pass from the model (1) to the spin-fermion model, 
which is justified in the vicinity of magnetically order state in 
two-dimensional systems J8fl. We employ the quasistatic ap- 
proximation for magnetic fluctuations in the effective field [81]. 
The static irreducible susceptibility reads 



d 3 S 



in|,(q|S signS z )+ 
+ ^n ± (q\S signS 2 ) 



exp(-iS^oo), (2) 



where Z^«j is a normalizing factor, 

1 .. f(e k -US)-f(e k+q + US) 

k 



n, u (q|S) 



-IlimV^ 

Nt^o^i e k - 



e k+q - US ± US 



S^oo = 31/ S /(2A ) is gaussian spin fluctuation field action, 
specified by only one parameter A (having a sense of dispersion 
of the uniform static mode S), which strongly simplifies the 
calculations. 



Minimal critical value U c , which is necessary for the stabil- 
ity of ferromagnetic ordering, is determined by a generalized 
Stoner criterion 

U e = l/// q=0 . (3) 
The r.h.s. of Eq. (0 can be simplified via the relation 



4+8f' 



(4) 



-4+8/' 



where p is bare DOS and ip(x) = (3x 2 + 2) exp(-3jc 2 /2)/ V67r. 

In fact, ferromagnetic state becomes stable when H q as a 
function of q acquires a maximum at the point q = 0: 



[H q=Q - H q ] > 0. 



(5) 



The application of the criterion (O together with Eq.© allows 
to calculate i/ c which determines the boundary of ferromag- 
netic and incommensurate magnetic phases in the ground state 
at arbitrary t' jt. 

Phase separation (PS) into Neel antiferromagnetic and ferro- 
magnetic phases may also play an important role in our prob- 
lem. Recently, the generalization of Visscher's method 11 1 Oil for 
PS boundary line has been performed for the case of square 
lattice Hi: 



tjU PS (n) = [1 + 2sign(l - n)t' jt]n(l 



■«) 2 /2, 



(6) 



One can see that PS boundary lines «ps(t/) possess strong 
asymmetry with respect to half-filling at finite t ' . 

In the limit of strong correlations we perform also calcula- 
tions of the electron Green's functions using the equation-of- 
motion method in the many-electron representation of Hub- 
bard's X-operators [11]. Note that within this approach we 
treat only the destruction of saturated or non-saturated ferro- 
magnetic ground state by one-particle (individual) spin-flip ex- 
citations. In the simplest Hubbard-I approximation the Green's 
functions read 



G ktT (E) 



SkY 



where the bare inverse locator F%.(E) has a form 



E(E - U) 
E - U(no + n<r)' 



(7) 



(8) 



n a being the number of holes (a = 0) or singly-occupied states 
with spin projection <x. However, this approximation is inap- 
plicable to the ferromagnetism problem since it violates impor- 
tant kinematic relations. It was improved in Ref. iflHl by in- 
cluding spin and charge fluctuation to obtain F® T (E) — > F ktr (E), 
where F kcr is the corrected inverse locator which is determined 
self-consistently. This many-electron (ME) approximation al- 
lows to calculate the quantities n a and obtain the equation 
for magnetization. Solving this equation yields two magnetic 
phase transitions: from saturated ferromagnetism (sFM) to 
non-saturated ferromagnetism (nsFM), and from non-saturated 
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Figure 1: Ground state magnetic phase diagram in variables (n,U). f = 0. 
Comparison of different types of instability of ferromagnetic ground state. "QS: 
FM-IC" is boundary line between ferromagnetic and incommensurate mag- 
netic phase regions obtained within quasistatic approach. "ME: sFM-uFM" 
is the boundary line between saturated and unsaturated ferromagnetic phase re- 
gions obtained within ME approach, "ME: uFM-PM" is boundary line between 
unsaturated ferromagnetic and paramagnetic phase regions obtained within ME 
approach flUl . "VS: PS region" is boundary region of PS to saturated ferro- and 
antiferromagnetic phases obtained within Visscher method (see Q|) 



ferromagnetism to paramagnetic state (PM). The first instabil- 
ity has in fact spin-polaron nature (occurrence of the pole of the 
spin-down Green's function below the Fermi level). 

The phase diagrams calculated by two methods described are 
shown at small t'/t in Figs. Q]and[2] We see that the spin- 
density wave (Eqs. (0 and (|5]l) and one-particle (Eqs. (0 and 
(JSJ) instability lines restrict the ferromagnetic region from dif- 
ferent sides. 

The ferromagnetic state for f — occurs for very large 
(practically unrealistic) U only because of competition with 
commensurate or incommensurate antiferromagnetism. With 
increasing t' , the van Hove singularity (VHS) in the electron 
spectrum, which determines ferromagnetic instability, is shifted 
from the band centre, and the critical U considerably decreases 
for n < 1 . Thus the next-nearest neighbor hopping is of cru- 
cial importance for ferromagnetism, the magnetic phase dia- 
gram being strongly asymmetric with respect to half-filling. 

3. Instabilities of ferromagnetic state: t'/t close to 1/2 

In the case of large ?' (in our calculations, ?' = 0.45f) ferro- 
magnetism is practically absent (at not too large U/ 1) for n > 1 . 
Hence we focus our attention on the case n < 1. When the 
Fermi level lies in the vicinity of VHS, two opposite tenden- 
cies take place: (i) the stability of ferromagnetic order substan- 
tially increases due to large DOS at the Fermi level (critical U 
decreases), (ii) long-wave spin-density instability (rather than 
the Neel antiferromagnetism) becomes actual ferromagnetism's 
competitor, as opposed to the case of small t'/t (critical U in- 
creases). 



Figure 2: The same as in Fig. f = 0.2? 



Since critical U necessary for the ferromagnetism stability is 
not large, we do not apply the strong-coupling ME approach 
used in previous Section. The results of the quasistatic ap- 
proximation |@] are shown in Fig. [3] A comparison with the 
results of Ref. [7] at t' w 0.43f is also shown. The phase 
diagram is strongly asymmetric with respect to VHS: when 
Fermi level lies above VHS, competition with long-wave mag- 
netic order increases critical U. At the same time, when Fermi 
level lies below VHS, ferromagnetic ordering is stabilized by 
rather small Coulomb interaction. The asymmetry is connected 
with the electronic topological transition which occurs when the 
Fermi level crosses VHS. At the Fermi level above VHS, flat 
parts of the Fermi surface occur in the vicinity of VHS points 
("quasinesting" situation), which results in enhancement of the 
incommensurate fluctuations. This conclusion is confirmed by 
the functional renormalization group studies [9]. We can also 
see that our results are in qualitative agreement with the results 
of the Gutzwiller-type method [7]. However, the latter does not 
capture delicate features of incommensurate fluctuations prop- 
erties, in particular, strong dependence on the Fermi level posi- 
tion. 

To conclude, we calculated the boundaries of ferromagnetic 
region on the ground state phase diagram at different t' . We 
found that ferromagnetic region is restricted by instabilities of 
two types: spin-wave instability and single-particle instabil- 
ity. It was shown that the Stoner criterion is not applicable for 
both small and large f . Indeed, at small t' in the vicinity of 
half-filling the ferromagnetism can be destroyed by collective 
spin-wave excitations with large wave vector. Away of half- 
filling the ferromagnetic order is destroyed by single-particle 
excitations, but such an (spin-polaron) excitation causes strong 
many-electron renormalizations and should be treated with ac- 
count of electron correlations. At large t'/t < 1/2 the unique 
possibility of ferromagnetic order occurs in the vicinity of van 
Hove singularity at relatively small U, but the long-wave spin- 
density instability substantially suppresses the ferromagnetic 
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Figure 3: Ground state magnetic phase diagram in variables {n, U). f = 0.45. 
Comparison of results obtained within quasistatic approach ("QS"), Stoner the- 
ory ("Stoner") and Gutzwiller-type approach ("Hanisch") | 7] is shown 



order in comparison with the Stoner criterion. 

The above results for the case of large f'/f yield a possi- 
bility to explain the magnetic properties of doped ruthenates 
like Sr2-.iLa t Ru04, where the ferromagnetic order appears to 
be suppressed [12] despite that Fermi level is close to VHS 
of flf rv -derived band and the Fermi energy is above VHS H 1 3fl - 
Since tight-binding fit for the latter band yields f'/f ~ 0.4 01411 
and x-ray studies estimate U as about 3.5? lfl5ll . we can ex- 
plain the ferromagnetism suppression by competition with in- 
commensurate long-wave magnetic order. Moreover, one can 
expect that hypothetical increase of available range of La con- 
centration will induce the transition into the ferromagnetic state 
at the x value where the Fermi energy coincides with VHS of 
(/^-derived band. Thus the difference in electron spectrum of 
layered systems (e.g., of cuprates and ruthenates) strongly in- 
fluences their magnetism. 
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